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Abstract 

We present an algorithm for the maximization of photonic bandgaps in two-dimensional 
crystals. Once the translational symmetries of the underlying structure have been imposed, our 
algorithm finds a global maximal (and complete, if one exists) bandgap. Additionally, we prove 
two remarkable results related to maximal bandgaps: the so-called 'maximum contrast' rule, 
and about the location in the Brillouin zone of band edges. 

1 Introduction 

Photonic crystals are structures where the electromagnetic properties of the constituent material are 
periodic in space. They have been intensely studied recently Q both for their practical applications 
and the theoretical challenges they pose. One area of research is the design and fabrication of 
structures that do not allow electromagnetic waves with frequencies in a certain range to propagate 
in them. Such structures are expected to find use as perfect mirrors and guiding structures, etc. In 
many applications, it is desirable that the excluded frequency range - the bandgap - be as large 
as possible. In this paper, we focus our attention on the design of structures with such maximal 
bandgaps. 

The design problem that we propose to solve has been addressed before j|, Hi @|< The starting 
point in each approach has been to impose the translational symmetries of the structure to be 
designed, and to impose global constraints that the material properties must staisfy (maximum and 
minimum dielectric constants, for instance). We do the same. In earlier approaches, the next step 
involved exploring a subset of the structures that satisfy the imposed translational symmetry. Our 
algorithm, in contrast, involves an unrestricted exploration of all the structures that satisfy the 
imposed conditions. Notably, the algorithm is polarization non-specific; it can maximize a complete 
bandgap (i.e. a frequency range that is excluded for all polarizations). 

The algorithm we propose is based on incremental steps, each of which increases the size of the 
selected bandgap. Each step is based on a variational argument that is discussed in Q and outlined 
below. 

The plan of the paper is as follows. In the next section, we present the equations relevant to 
the problem. Thereafter, we outline the variational argument we use and describe an algorithm to 
implement this in the case where the dielectric constant of the material is required to lie in a certain 
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range. We have numerically implemented the algorithm; the details of the implementation will be 
presented elsewhere. 

Two remarkable results follow quite easily from our discussion: First, that for a maximal bandgap 
structure where the dielectric constant is allowed to vary within a range, the dielectric constant at 
each point takes on either the smallest or the largest possible value. Second, that the wavevectors 
corresponding to the edges of a maximal bandgap may take on values only from a very small subset 
of the Brillouin zone. 



2 Maxwell equations in two-dimensional systems 

We consider a two-dimensional photonic crystal with a dielectric constant e(r) that is periodic 
in two independent directions (that define the x-y plane) and homogeneous in the third (the z- 
direction). The magnetic susceptibility of the constituent material is assumed spatially constant. In 
this system, electromagnetic waves with wave vectors in the x-y plane come in two distinct species: 
iS-polarization, where the electric field vector is parallel to the z-dircction; and, ii-polarization, 
where the magnetic field vector is parallel to the z-direction. In the former case, Maxwell equations 
yield 
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and in the latter, 
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E and H are the z-components of the electric field and magnetization vectors, respectively. (See B 
for details of derivation.) Each of these equations has the form Du — Aw, where D is a Hermitian 
operator. In the case of i?-polarization, 

and u = e 1 / 2 !?. For -polarization, 

D = (—e- 1 — + — e- 1 — } (4) 
\ dx dx dy dy J 

and u = H. In both cases, A = ~ui 2 /c 2 . For ease of exposition, we will set c = 1. 

It is well-known that the spectra of eigenvalues of the equations above display 'bands'. Given two 
bands, the difference between the smallest eigenvalue of the upper band and the largest eigenvalue of 
the lower band characterizes the bandgap. (We will focus our attention on the spectrum of u>, rather 
than A.) After having selected two adjacent bands, we will seek a function e(r) that maximizes the 
bandgap between them. In general, e(f) may be subject to imposed constraints. While the method 
outlined below may be used to handle a variety of constraints, we will confine ourselves to a specific 
one: e min < e(r) < e max . 



3 The variational algorithm 

Starting with the standard eigenvalue problem Du — Au, where D is a Hermitian operator, we will 
first investigate changes in A due to changes in D. Prefixing 5 to denote 'change in' 

(it + S u y(D + SD)(u + Su) = A + SX. 
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Retaining terms to first order in the changes, and noting that for normalized eigenfunctions, u'Su ■ 
0, 

u'SDu = SX. 

Change in D is due to change Se in e. In the case of .E-polarization, 

}2 pa 
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so that 
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J-lD-D 5 -^ (5) 
2e 2e V ; 

u'SDu = — X ( u* — udr, (6) 
Jc 6 

where we use dr to denote the volume element of the unit cell C with boundary S. 

One may proceed similarly in the case of i?-polarization. Using the notational shorthand V = 
id/dx + j d/dy, 

u^SDu = [ u*V • (5(e _1 )Vu) dr 
Jc 

(u*S(e- v )Vu) -hdA- J (Vu*) ■ (Vu)^(e^ 1 ) <ir 

^(e-^lVu^dr. (7) 



Recall that SX = —IloSlo. To sum up, we have the following: 

Suj ( ui\u\ 2 /2e : B — polarization 



Se \ — |Vu| 2 /2we : H — polarization. 

It is not w but, rather, a difference lo-i — u>±, that we seek to maximize with respect to variations 
in e. Thus, we would like to tailor variations in e so that S(u>2 — Wi) > 0. This is achieved if Se > 
wherever Su)2/Se > Su)%/Se, and vice versa. If e = e ma x (e = £mm)i it is impossible to achieve Se > 
(<5e < 0). Therefore, maximal — is expected when e = e maa; wherever Su>2/Se > Su)i/Se, and 
e = fmm wherever Su>2/ 'Se < Su)i/Se. Note that no explicit reference to polarization has been made; 
ui\ and ix>2 may even be eigenvalues corresponding to different polarizations. This makes it possible 
for our algorithm to maximize the size of complete bandgaps. 

This suggests the following iterative algorithm to a maximal bandgap: 

0. Choose periodicity of e(f), and the adjacent bands that straddle the bandgap to be maximized. 
Initialize with an arbitrary (but appropriately periodic) e(r). 

1. Determine, using e(r), the spectrum of eigenvalues to and the wavefunctions u. 

2. Determine the wavefunctions that correspond to the smallest eigenvalue of the upper band (u2 
and 0J2) and the largest eigenvalue of the lower band [u\ and u>\). 

3. Make small changes Se in e(r): Se > wherever Su>2/8e > Suji/Se (determined using the 
expressions in (@)), and vice versa. (Maintain consistency with the constraint e TO j n < e(r) < 

4. Go to 1 unless termination/convergence criteria are met. 
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A discussion of some of the finer points of each step follows. 

The initialization step requires the lattice parameters as input, both to establish a length scale 
and to fix the translational symmetry of the system. It is in this sense - and this sense only - that 
the maximum attained is not global. (Observe that, due to the scaling properties of the eigenvalue 
equations, maximization without an imposed length scale is trivial.) 

Considerable work Q has been done towards numerical solution of the eigenvalue equations 
for arbitrary e(r). The usual practice is to first discretize the space in some manner, then select a 
reasonably dense subset of wavevectors from the Brillouin zone and, finally, obtain the eigenvalues 
corresponding to each wavevector. 

Remarkably, it turns out that it is not necessary to perform a complete bandstructure compu- 
tation to implement our scheme. To see this, note first that, since e(r) is real, wavevectors ±k have 
degenerate spectra. Therefore, any arbitrary linear combination u = otu^ + /3u_£ is an eigenfunc- 
tion. Note, too, from our earlier discussion (equation (||) and following paragraph), that either of 
the wavefunctions u = u<x or u = u\ that straddle the bandgap must be such that \u\ 2 (or |Vit| 2 ) 
has the same translational symmetry as the underlying lattice. This severely limits the possible 
values of k at the band-edges. In the case of a square lattice, it is easy to see that the only possible 
values of k are those with or ±7r as components. It is, thus, possible to limit the search - for a 
maximal-bandgap producing configuration of e(r) in a 2D square lattice - to wavevectors 0, (0, 7r), 
(7r, 0) and (ir, ir)l 

Another remarkable result - referred to following equation (|J) and hinted at in step 3 of the 
algorithm - is the'maximum-contrast' rule: if e(f) is allowed to vary arbitrarily between two limits 
then maximal bandgaps are obtained when e(f) attains one or the other extremal value at each 
point. (See |i| for a related detailed discussion.) 

One question that we have not answered here is whether the iterative algorithm converges. 
In a one-dimensional system (where exact results may be obtained otherwise and where E- and H- 
polarization spectra are degenerate), the algorithm does indeed converge rapidly to expected results. 
In 2D, numerical results obtained so far provide very strong evidence in favor of convergence. These 
results will be presented elsewhere M . 
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